Charge-induced maximal spin states of a polynuclear transition-metal complex 
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We theoretically investigate the ground state spin of a polynuclear transition-metal complex as 
a function of the number of added electrons taking into account strong electron correlations. Our 
phenomenological model of the so-called [2 x 2] -grid molecule incorporates the relevant electronic 
degrees of freedom on the four transition-metal centers (either Fe 2+ or Co 2+ ) and the four organic 
bridging ligands. Extra electrons preferably occupy redox orbitals on the ligands. Magnetic inter- 
actions between these ligands are mediated by transition-metal ions and vice versa. Using both 
perturbation theory and exact diagonalization we find that for certain charge states the maximally 
attainable total spin (either Stat — 3/2 or Stat = 7/2) may actually be achieved. Due to the Nagaoka 
mechanism, all unpaired electron spins couple to a total maximal spin, including unpaired electron 
spins on the metal-ions in the case of Co 2+ . The parameters are chosen to be consistent with 
cyclovoltammetry experiments in which up to twelve redox states have been observed. The above 
effect may also be realized in other complexes with an appropriate connectivity between the redox 
sites. The maximal spin states of such a charge-switchable molecular magnet may be experimentally 
observed as spin-blockade effects on the electron tunneling in a three-terminal transport setup. 

PACS numbers: 



I. INTRODUCTION 



The synthesis and investigation of molecule based 
magnets have been active areas of research for al- 
most a decade Q, 0. Today molecule based mag- 
nets offer a variety of chemical and magnetic proper- 
ties which have potential applications in a wide range 
of systems. Photo-magnetic (e.g. K .4Coi.3[Fe(CN) 6 ]) 
or spin-crossover magnetic substances (e.g. Fe(o- 
phenanthroline)2(NCS)2) offer potential applications in 
switchable devices. Polynuclear transition-metal com- 
plexes are particularly attractive in the effort to de- 
sign molecules with magnetic centers which couple ei- 
ther ferro- or antiferromagnetically. Interesting quan- 
tum tunneling effects were reported for ferromagnetic 
molecules, e.g. Mni 2 3] and Fe§ 0. Antiferromag- 
nets are of interest also due to effects associated with 
the Neel- vector 0, 0, Q • One interesting class of highly- 
designable complexes is the \M x Ml-grid structure which 
is formed by self-assembly 0, |j| [H], ^J. It consists 
of M 2 transition-metal centers which are positioned by 
two perpendicular arrays of rod- like ligands, each with 
M coordination sites. The transition-metals typically 
lie in a plane due to symmetric positioning of the lig- 
ands. Several realizations of these supramolcculcs exist 
which exhibit a variety of magnetic and electrochemi- 
cal properties, which can be manipulated via sidegroup 
substitutions. In this paper we focus in particular on the 
[2 x 2] -grid depicted in Figure^ It consists of four organic 
ligands (bis(bipyridyl)-bipyrimidine) and four transition- 
metal ions (e.g. Mn 2 +, Fe 2+ , Co 2 +, Ni 2 + or Zn 2 + [T^). 
Each metal-ion is situated in an approximately octahe- 
dral environment of nitrogen atoms. The total charge 8+ 
of the complex is countered by anions, e.g. BFJ . The 
Co 2+ -[2 x 2] -grid exhibits an impressive twelve reversible 




FIG. 1: Structure of the [2 x 2]-grid-type complex. 



reduction steps in solution in a cyclic voltammetry at 
T — 253K R is also possible to trap and iden- 

tify these complexes on a graphite surface in a controlled 
way [l^. Furthermore, in Co 2+ -[2 x 2]-grids the spin of 
each metal- ion was found to be S = 1/2 at low tem- 
perature and their intramolecular magnetic coupling was 
demonstrated to be antiferromagnetic 0] . 
We have developed a phenomenological strongly- 
correlated electron model which represents the important 
spin and charge degrees of freedom that play a role in 
both the redox chemistry and the magnetic interactions 
of [2 x 2]-grid complexes. We include the strong elec- 
trostatic interactions between electrons localized on the 
metal-ions and the ligands as well as the weak tunnel- 
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ing of electrons between them. A crucial feature of our 
model is that the orbital and charging energies on the 
ligands and ions are very different, a situation common 
in transition-metal complexes. We are interested in how 
the electron spins on the molecule couple to form the to- 
tal magnetic moment and how this coupling is affected 
by total electron number which can be electrically con- 
trolled in a three terminal transport setup (which is not 
considered here). The metal-ions mediate the coupling 
between electrons on the ligands and vice versa. We also 
want to study the effect of an unpaired electron spin on a 
mediating metal-ion on the ligand-ligand spin coupling, 
which will in turn determine the total spin formed by the 
ion and ligand subsystems. Two simple cases are there- 
fore considered, assuming each ion is in a low-spin state: 
either 5 = 0, Fe 2+ , or S = 1/2, Co 2+ . The Coulomb 
interaction also generates direct exchange coupling be- 
tween unpaired electrons on ligands and metal-ions which 
we also take into account. We focus on how the spins 
couple to a magnetic moment of large magnitude due 
to strong Coulomb and kinetic effects depending on the 
total charge. We do not consider effects on smaller en- 
ergy scales which determine the preferred direction of the 
magnetic moment. 

As expected, an antiferromagnetic coupling between the 
metal-ion and ligand subsystems prevails for most to- 
tal charge numbers. Extra electrons added to the bridg- 
ing ligands introduce a ferromagnetic spin correlation be- 
tween the metal-ions if the latter have a spin. However, 
we find that near half-filling of the four ligand orbitals 
sufficiently strong Coulomb interactions can stabilize a 
maximal spin ground state (S to t = 3/2 for Fe 2+ and 7/2 
for Co 2+ ) by the Nagaoka mechanism [lq. Due to the 
Pauli principle a missing or excess electron on the lig- 
ands (relative to half-filling) can be delocalized maxi- 
mally when the "background" of the remaining electrons 
is completely spin-polarized. In our model this effect 
is most easily achieved for a Fe 2+ -[2 x 2]- grid. For a 
Co 2+ -[2 x 2]-grid the spins of electrons localized on the 
bridging metal-ions counteract the Nagaoka mechanism. 
However, the direct exchange coupling between unpaired 
electrons on the neighboring ligands and metal-ions co- 
operates with the Nagaoka mechanism 1 1 6| . The result- 
ing maximal spin is however more than twice as large as 
for Fe 2+ . The charge-sensitive spin-polarization requires 
a strong local Coulomb charging effect and metal-to- 
ligand charge-transfer (MLCT) barrier. The former can 
be achieved chemically by introducing electron-donating 
sidegroups on the ligands 01- The electrochemical ex- 
periments 01 were performed near room temperature 
where thermal occupation of high-spin excited states of 
the metal-ions becomes possible. Generalization of our 
model to incorporate the high-spin states of the ions is, 
however, nontrivial since vibrational degrees of freedom 
of the nuclear framework are involved in their stabiliza- 
tion 0. This is beyond the scope of the present paper 
where we focus on low temperature behavior. 
Experimental detection of sublattice and total magne- 



tization as a function of the number of added electrons 
would be of great interest. In particular, we propose elec- 
tron transport through the complex in a three terminal 
setup. In a gate voltage range where the ligand orbitals 
are near half-filling the large change in ground state spin 
(AStot = 3/2 (Fe 2+ ), resp. AS tot = 7/2 (Co 2+ ) ) may be 
revealed through the spin-blockade effect 18]. We have 
investigated this transport problem in detail for a sim- 
ple generic model n [lJj . The goal of the present paper 
is to derive and investigate the effective model, in par- 
ticular the effect of bridging metal-ions with a spin. Ex- 
periments in semiconductor heterostructures show that a 
transport current can couple through the spin-blockade 
effect to electron spins localized in a quantum dot and 
subsequently to nuclear spins in the local environment 
via the hyperfine interaction | 20| . This offers an interest- 
ing perspective for single-molecule transport where the 
immediate nuclear environment of the device can be con- 
trolled at the level of chemical synthesis. 
The paper is organized as follows. In section [H] we dis- 
cuss the minimal set of orbitals involved in electron addi- 
tion processes and the intramolecular magnetic coupling. 
This serves as a motivation for the correlated electron 
model which we introduce in section II I II In section IIVI 
we discuss a perturbative treatment of the ion-ligand tun- 
neling in this model to gain a simple understanding of the 
role of spin degrees of freedom on the ions. We obtain 
the ground states and low energy total-spin excitations in 
sections I V Al and IV Bl We discuss the electron addition 
energies and the modulation of the ground state spin as a 
function of the number of added electrons. We conclude 
with a discussion of implications for electron tunneling 
experiments in section IVll 



II. ELECTRONIC STRUCTURE 

We can expect that the dominant molecular orbitals 
(MOs) at the Fermi energy are either d-like orbitals 
from the metal 2+ ions or 7r-orbitals from the ligands. 
The electrochemical experiments jlJI indicated that ad- 
ditional electrons occupy orbitals localized mostly on the 
ligands. which is not uncommon for polypyridine com- 
plexes plf . 

Let us first consider the individual metal-ions which are 
coordinated by six nitrogen atoms. In a ligand-field pic- 
ture, an octahedral coordination of each metal-ion in 
the grid would cause their d-orbitals to split up into 
two shells, tig {d xy ,d xz ,d yz ) and e g (d x 2_ y 2,d z 2). How- 
ever, the symmetry of the local environment is lower and 
results in a small splitting of the levels in each shell. 
The different orbital occupations for Fe 2+ and Co 2+ are 
shown in Figure [21 In the case of Fe 2+ in its low-spin 
configuration the ti g shell is closed and the e g orbitals 
are empty. In the case of Co 2+ there is an additional 
unpaired electron in the lowest of the two e g orbitals. 
The highest occupied molecular orbitals of the ligands 
are mostly of 7r character as ab-initio calculations show. 
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FIG. 2: Orbital configuration of transition metal-ions and 
ligands. 
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FIG. 3: Geometry of transition-metal [2 x 2]-grid consisting of 
metal-ion orbitals (small circles), Fe 2+ (left) and Co 2+ (right), 
connected by ligands (large circles). The hopping t, orbital 
energies E, e and Coulomb charging U, V, u, v, w and direct 
exchange energy F are shown schematically. 



The important point is that the empty or singly occu- 
pied e g -orbitals on the metal-ions have different approx- 
imate symmetry (cr). This means we can consider lig- 
ands and metal-ions weakly coupled systems and discuss 
the tunneling of electrons between them. The remaining 
question is whether extra electrons prefer to occupy the 
metal-ion d-orbitals or the ligand 7r-orbitals. As far as 
charging effects are concerned, one expects that the on- 
site energy U on the metal ions is larger than the on-site 
energy u on the ligands because orbitals are more con- 
tracted on the metal-ions. Furthermore, electrons would 
favor positions between two positively charged ions, i.e. 
on the bridging ligand. This is in line with the interpreta- 
tion of the cyclovoltammetry experiments | 12| . The main 
difference in our model between the Co 2+ /Fe 2+ -[2 x 2]- 
grids will thus be the presence/absence of an localized, 
unpaired electron in an orbital of the metal-ion bridging 
the two ligands (Figures [21 0) ■ 



III. CHARGING MODEL 

We model the [2 x 2]-grid complex by four metal (Fe 2+ 
or Co 2+ ) and four ligand sites with one spin-degenerate 
orbital per site (Figure |3J) ■ The following Hamiltonian 
then captures the features of the electronic degrees of 



freedom discussed in section ITT1 

H = Hj + H\_ + Hf e /Q Q + Hdw + Hy, 

h t = Y,Y, tA l° a j>° +h - c - 

(id) ° 
4 

H L = ^{erij +un jA n jd +vn j n j+1 ) 

2 

4 

tfpe/co = ^2(ENi + UN iA N itl ) 



(1) 
(2) 



Hy 



Vj2N inj . 



(3) 
(4) 

(5) 
(6) 



Operators and variables (except t,F,V) in lower/upper 
case relate to the metal- ion/ligands and it is implicitly 
understood that all indices appearing run from 1 to 4 
(e.g. j + 1 — > 1 for j = 4). denotes a summation 

over nearest neighbor metal-ions (i = 1 — 4) and ligands 
(j = 1 — 4). The Fermion operator Oj a (a>j,cr) creates 
(destroys) an electron on ligand site j = 1 — 4 with spin 
projection a — ±1/2. The occupation number operator 
is defined as usual n^ a — aj CT aj )ff and rij = ^2 a nj, a . 
Similar definitions hold for the metal-ion (Ai }(T , iVj )(T = 

where r is the vector of Pauli matrices, is the electron 
spin of a metal- ion and Sj t k = | a , at a T a>IT 'Cbk,a' , is an 
operator related to the ligands. The tunneling term (0 
describes hopping between ligand and metal-ions. For a 
D2d symmetric molecular structure the hopping matrix 
elements t are independent of sites. In the following we 
choose t = 1 and all energies are expressed in units of t. 
The ligand-part of the Hamiltonian in © consists of a 
spin independent orbital energy e, the classical Coulomb 
repulsion terms on the ligand (u) and between adjacent 
(v) and the opposite ligands (w). Due to decreasing 
overlap with distance we have u > y/2 > w (e.g. for 
Fe 2+ u ~ 4:V ~ 0.3eV and v m 2w 12}). Equation gj) 
describes the isolated metal-ion orbitals with energy E. 
For the metal-ions we only consider the short-range in- 
teraction U because the orbital overlap between two ions 
is much smaller than that between two ligands. In or- 
der to describe the spin coupling of metal-ion and ligand 
electrons correctly, we include the direct exchange © 
with F > 0, which is known to stabilize ferromagnetic 
states even when it is weak ^(| (see sections llVI and Ivt . 
For consistency the metal-ligand charging energy |jfj| also 
needs to be incorporated. In general F < V are of the 
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FIG. 4: Energy diagram: Example of a charge transfer 
from one ligand to another ligand by virtual occupation of 
Fe 2+ /Co 2+ -ion state in the middle. Note that the MLCT 
barrier A = e — E is negative in the case of Fe 2+ and posi- 
tive in the case of Co 2+ . Also, for Fe 2+ the doubly occupied 
state with energy E + U can be neglected in the perturbation 
theory. 



same order. However, in the regime of interest the par- 
ticle number on the metal- ions is fixed (N — for Fe 2+ , 
N = 4 for Co 2+ , cf. section so the second term in 
(|5I6() yield a constant. In Figure |3] these interactions are 
schematically indicated. 

We study the parameter regime where the first eight ex- 
tra electrons will occupy the ligands as in the experi- 
ment ^3 • (To describe more than 8 reduction steps more 
than one orbital per ligand has to be taken into account. 
This is not our purpose here.) For Fe 2+ we must then 
assume that the two charge states of the ligand lie below 
the Fe 2+ orbital energy E: 



Fe 2+ : e < e + u < E < E + U 



(7) 



For Co 2+ the two charge states of the ligand lie between 
the singly and doubly occupied states of the metal-ion: 



Co 2+ : E < e < e + u < E + U 



(8) 



The orbital energy difference A = e — E is associated with 
metal-ion to ligand charge transfer (MLCT) between un- 
occupied metal-ion and ligand sites. In the case of Fe 2+ 
A < and in the case of Co 2+ A > 0. 



IV. PERTURBATION THEORY 
MODEL 



EFFECTIVE 



In the limit |A| > \t\ (Fe 2+ ), resp. A, U - A > \t\ 
(Co 2+ ) the charge transfer between ligands and metal- 
ions is suppressed. In order to gain qualitative insight the 
fluctuations of the orbital occupation around zero (Fe 2+ ) 
or one (Co 2+ ) can be treated using 2nd order Brillouin- 
Wigner perturbation theory (equivalent to a Schrieffer- 
Wolff transformation |22|). We thereby eliminate the 
charge degrees of freedom on the metal-ion sites and in- 
corporate their effect in an effective tunnel coupling be- 
tween the ligands. An example of a virtual process giving 
rise to this coupling is shown in Figure^ In the resulting 
effective model each metal-ion site is thus either entirely 



eliminated (Fe 2+ ) or characterized by a pure spin degree 
of freedom (Co 2+ ). This model contains the low energy 
properties of the mobile electrons on the [2 x 2]-grid. As 
a result we are left with the effective Hamiltonian (up to 
a constant): 

#Fe = EE T i^ + if L (9) 

(jk) ° 

< = E E {( K +J^H^ + jTt° s j° 

+ (K +(J- F)ar[)n^ + (J - F)t+°sJJ} 
+ H L 

= E E {( J - F5 3,k)TiS jtk + ^Ka] a a Ka } 



i j,k—i,i-\-l 



(10) 



The coupling constants T, K and J are 
2A 



T 



1 i 2 

2 V A U-A 
~ 2 ( A + U-A 



(11) 



For Fe 2+ the effective Hamiltonian is the extended Hub- 
bard model on four ligand sites with effective hopping 
matrix element T. In contrast, for Co 2+ we retain an 
eight site model: the effective Hamiltonian couples the 
spin and charge on the four ligands to the spin on the 
four metal-ions. The first two terms in I jlOj l describe 
tunneling between ligands with (J) and without spin- 
flip (K). The next two terms describe fluctuations of 
the charge and spin on the ligands. The exchange cou- 
pling J > 0, which is counteracted by the direct ferro- 
magnetic exchange F > 0, favors a correlated ground 
state where the ligand and metal-ion spins are coupled 
antiferromagnetically. This coupling always dominates 
over the tunneling amplitude K: J > \K\ ^ and K 
even vanishes for A = ^. The sign of K depends on 
whether A < %(K > 0) or A > %(K < 0) and deter- 
mines whether the amplitude K + J for the tunneling 
of electrons with spin a parallel to the local spin on the 
metal- ion Si is enhanced/suppressed relative to the am- 
plitude K — J for spin antiparallel to Si. 
We point out that the effective Hamiltonian (|10|l for 
Co 2+ contains no explicit interaction term between the 
metal-ion spins: all interactions are mediated by the 
electrons on the ligands. In order to describe the elec- 
tron addition effects on the metal-ion spin coupling our 
second-order perturbation theory suffices. In the absence 
of extra electrons on the ligands, only in fourth-order 
perturbation theory a weak effective antifcrromagnetic 
Heisenberg exchange interaction between the spins on 
the metal-ions appears. This superexchange is mediated 
by an empty intermediate ligand orbital. The four site 
Heisenberg-model with this effective coupling has been 
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studied in fljj and agrees with intramolecular coupling 
found experimentally. As soon as a ligand orbital con- 
tains one electron, the weak fourth order effect is super- 
seded by the second order coupling incorporated in the 
effective Hamiltonian (If Oil . 



V. ADDITION ENERGIES AND SPIN STATES 

We now present the results for the effective Hamilto- 
nian ll§)l. (|l Ofl (perturbation theory) and the full Hamil- 
tonian £Q) . We first study the addition energy spectra 
which reflect mainly the electrostatic effects and then fo- 
cus on the spin properties of the ground states and lowest 
lying excited states as a function of the number of added 
electrons n. 



A. Fe^+-grid 

To highlight the effect of the electrostatic interactions, 
we first consider the noninteracting limit of the effec- 
tive 4-orbital model © for Fe 2+ (u,v,w < \T\). Of the 
four eigenstates the lowest one lies at energy e — 2\T\, two 
states are orbitally degenerate at e (due to the 4-fold sym- 
metry axis) and the highest state lies at e+2|T|. Electron 
addition in this case would give rise to a pair of twofold 
degenerate peaks with a fourfold degenerate peak in be- 
tween. This is in clear qualitative disagreement with the 
experiments | l2j . The addition spectrum can only be un- 
derstood by including the interactions u > v > w > \T\ 
in our effective model J0J • The first electron reduces one 
of the four ligands. The next one goes onto the opposite 
ligand in order to minimize the Coulomb interaction. The 
third and fourth electrons reduce the adjacent ligands. 
For the next four electrons this sequence of processes is 
repeated, each time doubly occupying a ligand orbital. 
We thus have two sets of four reduction peaks separated 
by a large gap of order u. Each set of four consists of 
two pairs of closely spaced peaks (distance w) separated 
by a moderate gap 2v — w < u. The tunneling between 
equivalent ligands only weakly affects this picture. In 
Figure E] we plot the addition energies as a function of 
the tunneling amplitude t. 

Now we discuss the ground state spin as successive elec- 
trons are added to the ligands. Filling the levels in the 
noninteracting case of the effective model @ (u,v,w <C 
|T|) according to the Pauli principle the ground state 
spin is S to t = 1/2 for odd particle number n. For even 
n — 2,6 the ground state spin is S'tot = 0, whereas for 
half-filling (n = 4) S'tot = and S to t = 1 are degener- 
ate. In the presence of interactions charge fluctuations 
are suppressed. For u ^> \T\ this gives rise at n — 4 to a 
Heisenberg antifcrromagnet with a singlet ground state 
(Figure E). 

For sufficiently large u > (Figure 0) the ground state 
spin for odd n = 3, 5 is enhanced from the noninter- 
acting value Stot = 1/2 to the maximal possible value 



FIG. 5: Addition energy of Fe 2+ -[2 x 2]-grid as function of 
tunneling t (in units of w) for u — 15, v = 3, W = 1, A = 
-50, U = 100 calculated from the full model[T] For \t\ < |A| 
we are in the perturbative regime where no spin is localized 
on the Fe 2+ site and the effective model © applies. We re- 
produce the addition energy spectrum measured in 1121 where 
the spacings between the ground state energies with different 
n are due to electrostatic interactions. 




012345678n 

FIG. 6: Ground-state spin as function of the number of elec- 
trons n added to the ligands for u > Uth (dashed blue line 
= Fe 2+ , black line = Co 2+ without direct exchange) and for 
u < ttth (dotted red line in both cases). 



Stot = 3/2 (Figure EJ). Now the tunneling between the 
ligands plays a decisive role. Because double occupa- 
tion is suppressed, a single hole/electron (relative to the 
half-filled state n = 4) can maximally gain kinetic en- 
ergy when the background of the other electrons is fully 
spin polarized. This ferromagnetic alignment competes 
with the antiferromagnetic spin coupling due to superex- 
change processes. Which process dominates depends on 
the strength of the onsite repulsion u relative to the hop- 
ping \T\ i.e. oc \T\. The gap between maximal spin 
ground state and lowest excited state saturates at a value 
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0.05 - 




FIG. 7: Fe 2+ -[2 x 2]-grid: Ground-excited state gap as func- 
tion of u for n — 3, A = — 10, v = 2.25, w = 1. For 
u > u t h ~ 4.15 the ground state has maximal spin. 



~ 2\T\ independent of u due to the kinetic origin of the 
effect. This is the underlying mechanism for the Nagaoka 
theorem 0] which guarantees that the ground state has 
maximal spin if u is larger than some positive thresh- 
old value. It applies to the effective model because 
it fulfills a certain connectivity condition for the lattice, 
namely that a so-called "exch ang e loop" exists which is 
no longer than four sites [23, UM . This implies that all 
basis states with common S^ ot are connected with each 
other via nonvanishing matrix elements of @ 24|. In 
order to attain an observable effect one should have a 
moderate T on the one hand, and on the other hand 
the onsite interaction u must be enhanced with increas- 
ing T. The latter can be achieved by a chemical mod- 
ification of the ligands that draw charge density in the 
ligand LUMO orbitals. Taking typical parameters 0] 
|A| w leV,t w 10 -1 eV; A-E Nag w 10~ 2 eV, we estimate 
u±h ~ leV which is reasonable. 

We have checked that the interactions v, w increase the 
critical value for the Nagaoka state but do not de- 
stroy it 0] . We have also analyzed the effect of disorder 
by making the ligand sites inequivalent through different 
MLCT barriers Al. As expected, the Nagaoka state is 
stable if the change in the MLCT barrier A is smaller 
than |T|. Otherwise, due to the localization of electrons 
the Nagaoka effect, which is of kinetic nature, is sup- 
pressed. We expect disorder effects to be relatively weak 
since the ligands and the metal-ions form a highly sym- 
metric grid of equivalent centers. 



B. Co 2+ -grid 

The addition energy spectrum which we obtain for the 
Co 2+ -grid is qualitatively the same as for the Fe 2+ -grid 
since electrostatic interactions dominate. In the experi- 
ment the reduction peaks corresponding to the first 
three electrons are similar to those of Fe 2+ , in agreement 
with our model. However, the next five peaks exhibit 



roughly a constant spacing ~ 0.25 — 0.3eV, corresponding 
to the charging of one big "island" with better screening. 
This would require in our model to artificially change the 
parameters to u ~ v w w ~ O.SeV beyond n — 3. Obvi- 
ously effects become important which are not included in 
our electronic low temperature model, e.g. adding elec- 
trons could result in a change in the molecular geome- 
try which will lead to different electrostatic interactions. 
Also at the high experimental temperatures individual 
Co 2+ ions may be in the high-spin state where the tt- 
symmetric i2g-orbitals are singly occupied. These can 
couple more strongly to the ligand orbitals and increase 
metal-ligand charge transfer. Here we are interested in 
the low temperature regime however and assume low-spin 
(5 = 1/2) Co 2 + ions. 

The spin properties of the eigenstates of the effective 
model ff 1 01) without direct exchange interaction, i.e. F = 
0, are qualitatively similar to Fe 2+ (Figure |SJ) and will be 
considered first. We have a singlet ground state at half- 
filling (n = 4), and a Nagaoka maximal-spin ground state 
near half-filling (n = 3, 5) for sufficiently large charg- 
ing (U > u > uf h °). At half-filling n = 4 the antiferro- 
magnetic Neel-state has the largest weight in the ground 
state. The electron spins on the metal-ions couple ferro- 
magnetically due to the presence of electrons on the lig- 
ands and and vice versa. This is to be contrasted to the 
situation at n — where the metal-ion spins couple anti- 
ferromagnetically and at n — 4 for the Fe 2+ -[2 x 2]-grid, 
where the electron spins on adjacent ligands couple an- 
tiferromagnetically. The total spins of the metal-ion and 
ligand sublattice couple antiferromagnetically to a singlet 
ground state. The appearance of the Nagaoka state at 
n = 3, 5 has a different origin than in the Fe 2+ -[2 x 2]-grid 
since we have exchange loops longer than four sites. We 
do, however, have a bipartite lattice and hopping occurs 
only between the ligand and metal-ion sublattices. This 
is also a sufficient condition for the Nagaoka theorem to 
apply [Sllll. 

The four spins on the bridging Co 2+ metal-ions cause two 
quantitative differences from the case of Fe 2+ . Firstly, 
the maximal spin value attained in the Nagaoka state is 
simply larger, 5 to t = 7/2. Secondly, the threshold value 
of u^ for the appearance of this state is dramatically in- 
creased, u^ h ° « 10 4 in both the full and effective models. 
This is due to the antiferromagnetic exchange coupling J 
between the metal-ion and ligand sublattices. This sup- 
pression of ferromagnetism can however be understood 
by considering two extreme limits. In the limit where 
two sublattices are equivalent (A = 0, u = U) Nagaoka's 
mechanism is ineffective because the exchange paths are 
to long (> 4 sites, [24|). This implies that with decreasing 
A a higher M^ h ° is required to stabilize the Nagaoka state. 
In the opposite limit where the sublattices are well sepa- 
rated in energy (A, U — A ^> |i|), the effective model <|10[) 
applies and u^ h ° is decreased. The maximal spin state can 
still be achieved when U is sufficiently increased for fixed 
A, such that i < A < (7/2 we have K ea —J: tunneling 
of electrons with spins antiparallel to the metal spins is 
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- F=0.06 

- F=0.065 
-■ F=0.07 




FIG. 8: Co 2+ -[2 x 2]-grid: Ground-excited state gap as func- 
tion of u for n = 3, A = 10, U = 100, v = 2.25, w = 
1 and three different values for the direct exchange F = 
0.06, 0.065, 0.07. The non-monotonic behavior of the gap is 
due to several level crossings until the Nagaoka state is the 
ground state (beyond the rightmost zero). The direct ex- 
change "kicks" the system into the Nagaoka state: the thresh- 
old value is reduced when F — > J = 0.06. Even for smaller 
F > one still needs u > itf h ° to drive the system into the 
maximal spin ground state. 



suppressed, favoring a polarized ground state. However, 
the splitting AE is decreased due to the small effective 
coupling constants. 

The neglect of any (even small) direct exchange rep- 
resents an unbalanced treatment of the magnetic cou- 
pling, since the Nagaoka mechanism and the direct ex- 
change are known to cooperate to stabilize maximal spin 
states 0| . Even when it is weak, F w J, the direct 
exchange leads to a dramatic reduction of the threshold 
Coulomb energy for achieving the Nagaoka state to val- 
ues more close to those found above for Fe 2+ : ~ 5 
(Figure |SJ| and the gap between ground and excited state 
increases again. Furthermore it is possible that an en- 
hanced direct exchange stabilizes a maximal spin state 
even at half-filling n = 4. Still the spin can be switched 
from to 7/2 by going from n — 2 to n — 3 (or from 
n = 5 to n = 6). The excitation gap AE as function of u 
(Fig. [8} saturates at values of the order of 0.2 J. In this 
limit the direct exchange "kicks" ^(| the system into the 
Nagaoka state. For F ^> J the ferromagnetic coupling 
exceeds the exchange coupling and the resulting ground 
state is always ferromagnetic, independent of u indicating 
the Nagaoka mechanism is not relevant anymore. 



VI. CONCLUSION 

We have analyzed a strongly-correlated electron model 
for a [2 x 2]-grid complex with four transition-metal cen- 
ters (either Fe 2+ or Co 2+ ) to illustrate the interplay be- 
tween electron addition and intramolecular spin coupling. 
Our model contains both localized magnetic moments 
and delocalized electrons, in contrast to the customary 
description of molecular magnets. We have based our 
model on the addition energy spectra of the experiments 
in Ref.[l2^. the crucial input being that the extra elec- 
trons occupy ligand orbitals. We found large changes in 
the total spin, AStot > 1/2, of the molecule upon vari- 
ation of the total electron number due to the Nagaoka 
mechanism. The large charging energies on the ligands 
required for the high spin states can be tuned chemically 
by adding electron-donating groups to the ligands. Lo- 
calized spins on the mediating metal ions (Co 2+ ) coun- 
teract the Nagaoka effect, but the direct exchange cou- 
pling with the neighboring ligands can compensate for 
this. The total spin in the Nagaoka state for the Co 2+ 
complex is therefore more than twice as large as for the 
Fe 2+ complex. Low temperature electron tunneling ex- 
periments can access the change of the molecular spin 
as a function of added charge. Spin-blockade effects [Tg 
will dominate the single-electron tunneling around transi- 
tions between charge states with maximal spin |Tflj . Also, 
the iStot = 1/2 Kondo effect usually expected for odd 
n = 3, 5 will also be suppressed. In any case the Na- 
gaoka mechanism lowers maximal spin states in energy. 
Even as low lying total-spin excitations Nagaoka states 
have clear transport fingerprints due to spin-selection 
rules 0]. Other molecular complexes can also show the 
above behavior. It is essential that the connect ivity of 
the electron-accepting centers is appropriate |24| for 
the Nagaoka mechanism to be effective. 
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